Parametric Resonance and Cosmological Gravitational Waves 



Paulo M. SE 

Centra Multidisciplinar de Astrofisica - CENTRA and Departamento de Fisica, 
F.C.T., Universidade do Algarve, Campus Gambelas, 8005-139 Faro, Portugal 

Alfredo B. Hcnriquef0 

Centro Multidisciplinar de Astrofisica - CENTRA and Departamento de Fisica, 
Instituto Superior Tecnico, UTL, Av. Rovisco Pais, 1049-001 Lisboa, Portugal 
(Dated: December 17, 2007) 

We investigate the production of gravitational waves due to quantum fluctuations of the vacuum 
during the transition from the inflationary to the radiation-dominated eras of the universe, assuming 
this transition to be dominated by the phenomenon of parametric resonance. The energy spectrum 
of the gravitational waves is calculated using the method of continuous Bogoliubov coefflcients, 
which avoids the problem of overproduction of gravitons at large frequencies. We found, on the 
sole basis of the mechanism of quantum fluctuations, that the resonance field leaves no explicit and 
distinctive imprint on the gravitational-wave energy spectrum, apart an overall upward or downward 
translation. Therefore, the main features in the spectrum are due to the inflaton field, which leaves 
a characteristic imprint at frequencies of the order of MHz/GHz. 
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I. INTRODUCTION 

Although gravitational waves have not yet been di- 
rectly detected, those of cosmological origin are at 
present the object of an important research effort, as they 
will provide us with a unique telescope to the very early 
stages of the evolution of the universe. 

In this work we investigate the production of gravi- 
tational waves during the transition from the inflation- 
ary period to the radiation-dominated era of the uni- 
verse. We assume this transition to be dominated by 
parametric resonance, triggered by the coupling of the 
modes of an inhomogeneous bosonic field x to the homo- 
geneous inflaton field 0. We shall concentrate our efforts 
on those gravitational waves created by quantum fluctua- 
tions of the vacuum [J, |3, ■ This mechanism is indepen- 
dent from the generation of gravitational waves directly 
sourced by large inhomogeneities in matter distribution 
occurring during parametric resonance [i]. Although or- 
der of magnitude calculations suggested that the last 
mechanism would dominate the production of gravita- 
tional waves, during parametric resonance, we thought 
it would be of interest to settle the issue with an accu- 
rate calculation based on the quantum fluctuations of the 
vacuum. 

Parametric resonance may play a fundamental role in 
the reheating of the universe after the inflationary period 
0. Based on previous work on the subject we shall 
investigate a simplified model of parametric resonance 
which, nevertheless, keeps the most important features 
of the more complicated model. For instance, this is the 
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case when, as done in this paper, we neglect the backre- 
action from the metric perturbations and also perturba- 
tions of the inflaton field itself. This considerably simpli- 
fies the calculations, without missing the most important 
feature, namely, the sudden increase in the modes of the 
X field, due to their coupling to the coherent oscillations 
of the inflaton <j). 

We shall consider two different potentials for the infla- 
ton field, namely, a quadratic potential V{(j)) — m|0^/2 
and a quartic potential V{(l)) = A0^/4. For the coupling 
between the inflaton field (j> and the bosonic field x we 
shall use g^cfpy^ /2. The field x '^iH be taken to have 
a mass term V{x) — rn^x^ /2. We incorporate reheat- 
ing through an elementary decay mechanism of the two 
scalar fields into a relativistic radiation fluid character- 
ized by a density prad- The decay rates of the fields (j) 
and X Ei^re determined by the decay constants and F^, 
respectively. Arguments have been given that the x field 
is strongly suppressed before the transition period begins 
. The enormous increase of the scale factor during the 
inflationary era results in an exceedingly small value for 
X at the beginning of the transition. We are then justified 
in treating the field x &s arising from its own quantum 
fiuctuations, without any primordial classical part. 

Gravitational wave production will be calculated by 
the use of the Bogoliubov coefficients as continuous func- 
tions of time, which obey differential equations first intro- 
duced by Parker [l| (see also Ref. [Q] for a different deriva- 
tion of these differential equations). Besides being a the- 
oretically correct method of calculation of the graviton 
production, it has advantages over the frequently used 
sudden transition approximation. Associated with the 
sudden transition there is always an overproduction of 
gravitons of large frequencies, requiring the introduction 
of an explicit cut-off for frequencies above the rate of ex- 
pansion of the universe. This problem is automatically 
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solved by the use of the continuous Bogohubov coeffi- 
cients [9|. During the inflationary period, the potential 
term of the inflaton field dominates over the kinetic term, 
yielding a situation that closely resembles the behavior of 
a fluid with an equation of state of the form p = — p, also 
characterizing a cosmological constant. This, plus the 
fact that gravitational waves produced during the infla- 
tionary period will be completely red-shifted, simplifies 
the calculations, making it reasonable to begin our nu- 
merical simulations at the end of the inflationary period, 
with simple initial conditions for the Bogohubov coeffi- 
cients appropriate to a cosmological constant. 

This paper is organized as follows. In Section [IT] we 
describe the method used to calculate the gravitational- 
wave energy spectrum, present the equations to de- 
termine the continuous Bogohubov coefficients and de- 
rive appropriate initial conditions for these equations. 
The simplified model of parametric resonance, describ- 
ing the transition between the inflationary period and 
the radiation-dominated era, as well as the equations de- 
scribing the evolution of the universe from the end of the 
transition period till the present time, are introduced in 
Section [IIII Our numerical simulations, for two different 
potentials of the infiaton field (quadratic and quartic), 
are described in Section IIVI Finally, in Section |Vl we 
present our conclusions. 



II. COSMOLOGICAL GRAVITATIONAL WAVES 

We assume a flat universe and write the perturbed met- 
ric, in conformal time 77, as 

ds'^ = a^{ri) {-drf' + -|- /iy(77,x)] dx'-dx^] , (1) 

where the tensor perturbations hij are expanded in terms 
of plane waves 



(2^)3/2a(r7)V2fc 

X [ap(77,k)e,,(k,p)e*-^$(?7,fc) +H.C.] . (2) 

In the previous expression, G is the gravitational con- 
stant, p runs over the two polarizations of the gravita- 
tional waves, fc = |k| = l-najX = auj is the co-moving 
wave number, Op is the annihilation operator, £ij is the 
polarization tensor, and the mode function ^ obeys the 
equation of a parametric oscillator. 



+ ( fc - - ) ^ = 0, 



(3) 



where a prime denotes a derivative with respect to con- 
formal time. Under appropriate circumstances, we may 
have exponentially growing solutions, with the gravita- 
tional fleld pumping energy into the gravitational waves. 
Quantum mechanically this corresponds to graviton pro- 
duction. We shall calculate the amount of gravitons pro- 
duced using Parker's method 



In a expanding universe the ground-state and the op- 
erators change with time. 
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meaning that the annihilation ap(r\,V) and creation 
a|,(?7, k) operators also change, such change being codified 
in terms of time-fixed annihilation Ap(k) and creation 
v4^(k) operators through the Bogohubov transformation 



ap(7?, k) = fc)Ap(k) + /3*(77, fc)4(k). 



(5) 



In the previous expression, the Bogohubov coefficients a 
and (i obey the condition jap — = 1. Comparing 
any matter field expansion at two different times, using 
Eq. (O and the corresponding one for 0^(77, k), it can be 
shown that the coefficient fi gives the number of gravitons 
created, = {Nk{v))- 

To find the power spectrum, P{uj), we use the defini- 
tion of the density of states, uj'^duj/ (2Tr'^c^), and the fact 
that each graviton contributes with two polarizations, 
2hLL!; then, from the definition of the energy density in 
terms of the power spectrum, dE ~ P{Lu)duj, we obtain 



Pico) 



(6) 



Taking into account that the gravitational-wave energy 
density is given by pew — J we finally get the 

dimensionless relative logarithmic energy spectrum of the 
gravitational waves at the present time 779: 



1 



dpGW 



Pcrit(»yo) dlnuj 



(7) 



where pcrit is the critical density of the universe. 

As we have just seen, in order to obtain the 
gravitational-wave spectrum we need to compute the Bo- 
gohubov coefficient /3. This is done by solving the set of 
differential equations P, ^ 



a 



13' 



2k 



2k 



/3 + ae~2''=(''"''') 



which, upon the redefinition 

a = i(X + y)e*'^'(''-'''), 



become 



X' 



-ikY, 



Y' = -1 fc2 _ ^ 
fc V a 



X. 



(8) 
(9) 

(10) 
(11) 

(12) 
(13) 
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Because we are interested in graviton production dur- 
ing the transition between the inflationary period and 
the radiation-dominated era, the above set of differential 
equations should be integrated from the end of the infla- 
tionary period right up to the present time. This requires 
us to specify a model of evolution of the universe for the 
entire period under consideration in order to determine 
a!' I a (see next section) and also to specify, at the end 
of the inflationary period, the initial conditions for X(r\) 
and Yiji). 

Clearly, the set of differential equations p2)l and (fT3|) 
should be solved numerically. However, we can take ad- 
vantage of the fact that it admits an exact analytical solu- 
tion for the case of a de Sitter universe, in order to specify 
appropriate initial conditions for ^(77) and Y (r\). Indeed, 
during the inflationary period, the potential term of the 
inflaton field dominates over the kinetic term, yielding 
a situation that closely resembles the behavior of a cos- 
mological constant. Since, in this case, the scale factor 
is given by a(r\) = [H{rii — rj)]^^, where 77 < 771 (iji is 
an arbitrary constant and H is the Hubble parameter), 
Eqs. (ini) and ((13|) yield the solution 



+ C2 



Hm - v) 



JHvi-v) 



+ C2 



Cl 



-1 



1 



i 1 



(14) 

^ifc(jji-r)) 



-ifc(?7i-?7) 



(15) 



The condition \a\^-\(}\^ = 1 imposes a constraint on the 
integration constants ci and C2 , namely, c\ — c\ = 1 . Let 
us choose Cl = 1 and C2 = 0. Then, at the end of the 
inflationary period, 



X 



1 + 



k 



fc2 



.(17) 



These expressions will now be used as initial conditions. 

Note that a slow-rolling inflaton field, as the one in 
our model (see next section), does not mimic exactly a 
cosmological constant and, therefore, the Hubble param- 
eter is not constant during inflation, as assumed above. 
However, during the slow-roll phase, the Hubble param- 
eter decreases so slowly, that Eqs. ([TS]) and p7)) . with 
H = H{rii), constitute a very good approximation. 

To finish this section, let us point out that measure- 
ments of the cosmic microwave background radiation 
impose an upper limit on the energy spectrum for an- 
gular frequencies corresponding to the present size of 
the horizon, /iQriGvv(a;hor, ?7o) ^ ^ 10~^^, where Hq = 
F(2a}/(100 km Mpc"^s-i) andwhor = 2xl0-"/io rad/s 
. Timing observations of milliseconds pulsars and 
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Doppler tracking of the Cassini spacecraft also provide 
bounds on the gravitational-wave energy spectrum, re- 
spectively, /i§riGw(ii^ui, ?7o) < 2.0 X 10~® for Wpui = 
2.5 X 10"** rad/s ^ and /i§f^Gw(t^Cas, %) < 0.014 
for wcas = 7.5 X 10^^ rad/s 121. Foi' higher angu- 
lar frequencies, of the order of a few hundred rad/s, 
an upper limit is available from the Laser Interferom- 
eter Gravitational Wave Observatory (LIGO), namely, 
hl^Gw{^^,'no) < 3.4 X 10~^ [l3|. Finally, an integral 
bound can be derived from standard Big Bang nucle- 
osynthesis, Hq ilG-w{^^,'ilo)duj/ijj < 5.6 X 10^^, where 

ujn « 10"^ rad/s"!, [3. 



III. PARAMETRIC RESONANCE 

Assuming that the transition between the inflationary 
and the radiation eras is dominated by the phenomenon 
of parametric resonance, the Einstein equations for the 
transition period are written in the form 



BttG 



-a ptot, 



a 
a 



AirG 



a?{ptot - 3ptot), 



(18) 
(19) 



where the total energy density and pressure are given by 



1 

Ptot — 



Prad, 



Ptot — 



^ 2 2/ 2\ /2 

-a g ix )4> 



rPrad, 



(20) 



(21) 



where V{4>) represents the quadratic and quartic poten- 
tials, and g"^ {x^) (jP' / 2 the coupling between the two fields 
(j) and X- 

The equations for the inflaton and the rclativistic ra- 
diation fluid are, respectively. 



-aV^(p' (22) 



and 



(23) 



Note that when terms involving x^ appear in homoge- 
neous equations, like the ones above, the ensemble aver- 
age 
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(27r)3 



(24) 



should be used [14| . The same applies to terms involving 
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The non-homogeneous equation for the field x is writ- 
ten in the K-component form, for the mode functions Xk 



+ 2-x:, + {n^ + a'ml + a'g^cf') x. - -aV^^^. (25) 



Finally, the above set of equations is complemented 
with the equations describing the evolution of the uni- 
verse from the end of the transition period till the present 
time, namely. 



a 
a 



A-kG 



Pmat(?yo) 



+ Pmat(?7o) 



-I- (3W + l)/9do(7?0 

SttG 2 / \ 

Prad(%) 



a(??o) 



0(770 



3(1— ii))' 



,(26) 



3{l-w)' 



,(27) 



where for the dark energy we take w = 0.78 and the den- 
sity of radiation, matter and dark energy at the present 
time are, respectively, Prad(77o) = 4.6 x 10"'^^ kg/m'^, 
Pmatim) = 2.6 X 10^2^ kg/m^ and pddm) = 6.9 x 
10-27 kg/m^ 

The above set of equations (fT8 |) -([27 l) . together with 
Eqs. and for the Bogoliubov coefficients, are 
solved numerically in order to obtain the energy spectrum 
of the gravitational waves. 



IV. NUMERICAL SIMULATIONS 

Our numerical simulations begin at the end of the in- 
flationary period and continue up to the present epoch. 
We use a Runge-Kutte method with variable step to solve 
the system of differential equations (fTO|) - (P51) . Equations 
(HH]) and ([77]) are used to check the accuracy of the nu- 
merical solution. The integral in Eq. p4|) is computed 
using the Simpson rule, with a cut-off at k = 2'KaH and 
a interval of integration divided in 20 segments of equal 
length h = (7r/10)a'/a. Having determined the time evo- 
lution of a" /a, we then solve numerically Eqs. and 
(again with a Runge-Kutte method with variable 
step) for different values of a; = k/a, compute /S^ with 
Eq. pT|) and, finally, obtain nGw{^,Vo) from Eq. ([7]). 

The energy spectrum of the gravitational waves, 
^Gw{^,Vo)j is computed for values of u in the interval 
'j-'min < w < Wmax- The lowcr limit corresponds to a gravi- 
tational wave with a wavelength equal, today, to the Hub- 
ble distance: uj^^in = 27rc/dHubbio('7o) ~ 27r77(77o) = 1.4 x 
IQ--'^^ rad/s, where the Hubble parameter at the present 
time was taken to be H{rjQ) = 71 km s-^Mpc"^. The up- 
per limit corresponds to a gravitational wave which, at 



the beginning of the transition period between the infia- 
tionary and radiation-dominated eras, had a wavelength 
equal to the Hubble distance at that time; because of 
red-shifting, the angular frequency of this gravitational 
wave is, today, Wmax ~ '2TrH{ri^)[a{r]i)/a{riQ)], where 77^ is 
the conformal time at the beginning of the transition pe- 
riod. Within the model we are considering, the maximum 
angular frequency is about (10^ — 10^°) rad/s. 

There are two more characteristic values of uj in the 
gravitational-wave energy spectrum, corresponding to 
gravitational waves which had a wavelength equal to the 
Hubble distance at the moments 770->rad and Tyrad^mat, 
these being, respectively, the time when the energy den- 
sity of the inflaton becomes equal to the energy density 
of radiation (marking the end of the transition period be- 
tween the inflationary and radiation-dominated eras) and 
the time when the energy density of radiation becomes 
equal to the energy density of matter (corresponding to 
the transition from a radiation-dominated to a matter- 
dominated universe). These angular frequencies are, to- 
day, w^^rad « 10'' rad/s and Wrad^mat ~ 10"^^ rad/s. 

The gravitational-wave energy spectrum is then natu- 
rally divided in three regions. The first, from Wmin to 
'^rad^mat, marked by a sudden decrease of ^gw, the 
second, from ixirad^mat to [x)0_ti-ad, whcrc r^GW is con- 
stant (corresponding to a radiation-dominated universe, 
in which there is no production of gravitational waves), 
and the third, from c^J^^rad to cij,nax, where f^GW has some 
complex structure (peaks and valleys) , which depends on 
the details of the transition from the inflationary period 
to the radiation-dominated era. 

Our investigation is carried out for two different poten- 
tials of the infiaton field cj), namely, V{(j>) = m?^(fP' /2 and 
V[cf)) — A(/)'*/4. The values of the parameters and A 
are chosen such that the bounds from the temperature 
anisotropy of the cosmic microwave background radia- 
tion (CMBR) and from large-scale structure (LSS) are 
satisfied ^15j . 



and 



"^^ = 2iV(fc.) + l "^P' 



37r2p,(fc, 



(28) 



(29) 



In the previous expressions, Ps{k) is the power spectrum 
for density perturbations and N{k) is the number of e- 
foldings of expansion between the time when a co- moving 
distance scale labelled by k exited the horizon during in- 
flation and the end of inflation, both Ps(k) and N{k) be- 
ing evaluated at the CMBR/LSS scale k^ = 0.05 Mpc"^ 
For 47 < iV(fcc) < 62 and Ps{kc) = (2.45 ± 0.23) x 10"^ 
[Til l, we obtain the constraints 1.2 x 10~^ < m^/nipi < 
1.7 X 10-s and 1.3 x IQ-^^ < A < 3.6 x 10~^^. 

In the slow-roll approximation, inflation ends when the 
parameter e = [m'^i/ {16n)]{d^V/V)^ ~ 1, implying that 
0('7i) ~ mpi/{2y/n) for the case 1^(0) = rn^^cfp /2 and 
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FIG. 1: (Color online) The jump in the value of p^. Af- 
ter a short time interval, the energy density of the field x 
(lower/blue curve) becomes comparable to the energy density 
of the inflaton (upper /red curve). Both curves refer to case 
with 3 = 6.5 X 10"^. 
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FIG. 2: (Color online) Gravitational-wave spectra for the 
model with V^(0) — 771^0^/2. The effect of switching on the 
resonant field x is an overall translation upwards or down- 
wards. Upper/red curve refers to the case with g = 6.5x 10~^, 
middle/blue curve to (7 = (field % switched off) and 
lower/green curve to g = 10~^. 



^{r\i) ~ mpi/^/TT for the case V{(f>) — A0'*/4. The ini- 
tial value of 0' is chosen such that 0'^(r/i)/[2a^(?7i)] < 
F[(/)(?7i)]. We also choose XniVi) and x'niVi)! '^ell as 
and g, such that ptot [see Eq. (j20l) ] is dominated by 
the potential and kinetic terms of the inflaton fleld. Since 
any pre-existing radiation fluid is dilute during inflation, 
we choose PradiVi) — 0- Finally, the decay constants 
and are chosen to be small enough to allow a signifi- 
cant growth of the resonant field x before the evolution 
of the universe becomes dominated by radiation. 




A. Case V{(l>) 



1 1.5 2 
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The masses of the fields (j) and Xi the coupling con- 
stant g and the decay constants T^ and are cho- 
sen to take the following values: = 1.5 x 10~^mpi, 
= 10~^ mpi, 0<g< 10"^ = Tx = 5 X IQ-^" mpi. 
As initial conditions, at the beginning of the transi- 
tion period, we choose 0(77^) — 1, (t>{rii) = O.Srripi, 
(b'im) = -4.5 X lO-^m^j, p{t^,) = 0, Xn{m) - 10"'' 
and x'kI'?!) — the constraint equation (jlSp gives then 
a'{Tj^) « 1.30 X 10-6 mpi. 

The relevant parameter to consider is the resonance 
parameter, defined as g = g'^ (j)^ / {Aw?^ ^g), with the effec- 
tive mass of the inflaton given by = -I- g'^ix^)- 
In order to achieve a significant resonance we need values 
of 9 > 1, the so-called broad resonance regime. 

Our numerical simulations show that the energy den- 
sity of the field X, = (x'^)/(2a-^)+TO^(x^)/2, increases 
by many orders of magnitude in a short time interval, 
reaching values comparable to the values of the energy 
density of the inflaton field, = (j)'"^ / {2a?)+vn?^4>'^ /2. In 
Fig. 1 we can see that jump in the value of p^^. However, 
when switching off the resonant field x^ the most rele- 
vant change occurring in the spectrum f^cw is an overall 
translation upwards or downwards (see Fig. 2). 



FIG. 3: The change in time of the effective mass of the infla- 
ton. Its value increases sharply, when the energy density of 
the resonant field x becomes comparable to the energy density 
of the field (j). 



No explicit and distinctive signal left by the resonance 
can be read from f^GW- At first we assumed that this 
was due to the fact that p^^ did not dominate over p^ 
and prad for a sufficiently long period of time, and that, 
by increasing g (and, therefore, the resonant parameter 
q), this could be solved. This does not happen and, from 
Figs. 3 and 4, we may understand why: when (x^) begins 
to increase, the same happens to the effective mass of 
the infiaton, due to the presence of the term g'^{x^); as a 
consequence, the value of q decreases and the resonance 
is killed. 

As expected, the transition between the inflationary 
period and the radiation-dominated era leaves its imprint 
in the gravitational-wave spectrum ficw at frequencies 
of the order of MHz/ GHz. However, this imprint is due 
to the form of the potential of the inflaton fleld and to 
the values of the parameters and T^; the presence 
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FIG. 4: The change in the resonance parameter q. Because 
of the sharp increase of the effective mass of the inflaton, 
the resonance parameter suddenly decreases to values much 
smaller than the unity, q <C 1; as a consequence, the resonance 
regime ends. 
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FIG. 6: (Color online) Gravitational-wave spectra for the 
model with V{cj)) ~ A<^*/4. Upper/red curve refers to the 
case with g — 10"^ and lower/blue curve to g — (field x 
switched off). 
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gxlO^ 

FIG. 5: A manifestation of the bands of resonance in the 
model with V{(l>) = m%(j? /2. Q. Gw,o denotes the energy spec- 
trum of the gravitational waves in the absence of the resonant 
field X- The ratio f2Gw/f2Gw,o is taken at the flat part of the 
spectrum, namely, at tj = 10^ rad/s. For small values of g 
the spectrum tends to move downwards relatively to the case 
in which the field x is switched off, while for large values of g 
we have the opposite effect. 

of the resonant field x only induces an overall upward 
or downward translation of the gravitational-wave spec- 
trum. This upward or downward translation depends on 
the values of the mass of the resonant field , the cou- 
pling constant g and the decay constant T^ . In Fig. 5 we 
compare the energy spectrum in the presence of the reso- 
nant field X (for < <? < 10~^) with the energy spectrum 
in the absence of this field. The comparison is made at 
the flat part of the spectrum, namely, at w = 10^ rad/s. 
As clearly seen in Fig. 5, the spectrum tends to move 
downwards for small values of g, while for large values of 
g it tends to move upwards. 

The results presented above indicate that, as far as the 
creation of gravitons due to quantum fluctuations of the 
vacuum is concerned, the main features in the spectrum 



tjo) are due to the inflaton fleld. This is different 
from what happens when we consider gravitational waves 
directly produced by the motion of matter. In that case, 
the time-dependent inhomogeneities in the distribution 
of matter, produced by the rapid growth of the resonant 
fleld X, leave a much stronger imprint in the energy spec- 
trum of the gravitational waves flew at frequencies of the 
order of MHz/GHz 

B. Case V{<p) = 

The parameters are chosen to take the following val- 
ues: A = 3.5 X 10-"^ = rupi, < 5 < 2.8 x 10-^ 
= = 3 X 10~^^. As initial conditions, at the be- 
ginning of the transition period, we choose 0(77^) = 1, 
(t){rji) = 0.6 mpi, (t)'{r]i) = -1.5 x 10""^ m^^, p{rji) = 0, 

XniVi) = 10~^^ and XniVi) — 0; the constraint equation 
(froi gives then a'{r]i) « 4.35 x 10"'^r7ipi. 

In Fig. 6 we show two spectra, one with and the other 
without the field x- The same conclusions can be drawn 
as presented in the previous subsection: no obvious im- 
print of the resonance of the field x in the energy spec- 
trum of the gravitational waves, apart from the overall 
displacement. 

In Fig. 7 we show the bands of resonance which, in this 
case, are much better marked due to the fact that q is 
now the ratio between constant parameters, q = g'^ /{AX). 

If we compare Figs. 2 and 6, something interesting can 
be observed. While in Fig. 6 we see the oscillations in 
the MHz/GHz region to be around the value defined by 
the flat part of the spectrum, in Fig. 2 there is flrst 
a pronounced decrease, followed then by the usual os- 
cillations. This can be qualitatively understood as fol- 
lows. Given an homogeneous scalar field with a potential 
V{(f>) oc (/)", oscillating rapidly relatively to the expan- 
sion rate of the universe, for n = 2 the energy density of 
the scalar-field oscillations behaves like non-relativistic 
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FIG. 7: The structure of bands of resonance in the model 
with V"(0) = A<^''/4. n Gw,o denotes the energy spectrum of 
the gravitational waves in the absence of the resonant field 
X- The ratio Ogw/^^gw,o is taken at the flat part of the 
spectrum, namely, at lu — 10^ rad/s. 

matter (with p = 0), while for n = 4 it behaves like rela- 
tivistic matter with equation of state p = p/3 [l6l |. Now, 
it is known that the first case leads to a decrease in 
the energy density of the gravitational waves, while the 
other tends to keep the energy density at the level of the 
horizontal part of spectrum. 

V. CONCLUSIONS 

In this work we have investigated the production of 
gravitational waves due to quantum fluctuations of the 
vacuum during the transition between the inflationary 
period and the radiation-dominated era of the universe, 
assuming that the phenomenon of parametric resonance 
dominated this transition. The energy spectrum of the 
gravitational waves was calculated using the continuous 
Bogoliubov coefficients. This approach automatically 
solves the problem of overproduction of gravitons at high 
frequencies which afflicts the sudden-transition approxi- 
mation. 

We have found that the resonance field x leaves no ex- 
plicit and distinctive imprint on the gravitational-wave 
energy spectrum, apart an overall upward or downward 
translation. The amplitude of this translation, which de- 
pends on the values of the mass of the resonant field m^, 



the coupling constant g and the decay constant F^, is 
quite modest. Consequently, as far as the production of 
gravitons due to quantum fluctuations of the vacuum is 
concerned, the main features in the spectrum are due 
to the inflaton field, which leaves a characteristic im- 
print at frequencies of the order of MHz/GHz. For the 
models of chaotic inflation we have considered (quadratic 
and quartic potentials) , the relative energy density of the 
gravitational waves in this frequency range is quite small, 
r^GW ^ 10~^^. Therefore, our calculations show that, by 
far, the main contribution to the production of gravitons 
in the MHz/ GHz frequency range, for this type of mod- 
els, does come from the direct coupling of the anisotropic 
stress tensor, describing the motion of particles during 
preheating, as investigated in Ref. [3j. 

The MHz/ GHz region of the spectrum lies far beyond 
the range of frequencies accessible to ground- or space- 
based interferometer detectors, as well as to resonant bars 
and spheres (spanning roughly from 10~^ to 10"* Hz). In 
recent years, considerable effort has been made in or- 
der to develop high-frequency gravitational wave detec- 
tors [Tsj . In a near future, these detectors may reach 
sensitivities high enough to explore this region of the 
gravitational-wave spectrum. 

Below frequencies of the order of MHz, the influence of 
the transition regime between the inflationary period and 
the radiation-dominated era should not be felt anymore 
and, unless strong anisotropics develop later on, the main 
source of gravitational waves will then be the quantum 
fluctuations of the vacuum. However, the level of the flat 
part of the spectrum is mainly flxed by the parameters 
defining the kinetic and potential terms of the infiaton 
field, which are known to obey severe constraints from 
the observations. What our calculation indicate is that 
the two popular models of inflation we have used will not 
be able to generate enough gravitational waves to be seen 
by the interferometer detectors LIGO, Virgo and LISA 
or in the cosmic microwave background radiation. 
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